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Abstract 

The equation of state of hot hadronic matter is obtained, by taking into 
account the contribution of the massive states with the help of the resonance 
spectrum r(m) ~ m 3 justified by the authors in previous papers. This equation 
of state is in agreement with that provided by the low-temperature expansion 
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for the pion intracting gas. It is shown that in this picture the deconfinement 
phase transition is absent, in agreement with lattice gauge calculations which 
show the only phase transition of chiral symmetry restoration. The latter is 
modelled with the help of the restriction of the number of the effective de- 
grees of freedom in the hadron phase to that of the microscopic degrees of 
freedom in the quark-gluon phase, through the corresponding truncation of the 
hadronic resonance spectrum, and the decrease of the effective hadron masses 
with temperature, predicted by Brown and Rho. The results are in agreement 
with lattice gauge data and show a smooth crossover in the thermodynamic 
variables in a temperature range ~ 50 MeV. 

Key words: hot hadronic matter, resonance spectrum, equation of state, quark-gluon 
plasma, chiral symmetry restoration 
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1 Introduction 

One of the main goals of experiments with high-energy nuclear collisions is to produce 
and to study hadronic matter, in particular, trying to reach conditions at which the 
phase transitions into the quark-gluon plasma phase can take place [1-5]. The physics 
of hot hadronic matter has not been studied much, although such matter is already 
produced in current experiments in high-energy physics. 

The experimental data on multiple hadron production obtained in recent years are 
in agreement with the main consequences of the theory formulated by Landau || over 
40 years ago. However, with a sufficiently quantitative approach, it becomes necessary 
to consider a number of physical effects which bring about certain modifications 
of the results obtained in the fundamental work [6]. For instance, the solution of 
the equations of motion obtained by Landau differs from a numerical calculation of 
Milekhin 0, as a result of an inaccurate estimate. A more accurate analytic solution 
was given in |§. 

The equation of state in Landau's work is taken to be p = p/3 (where p is the pres- 
sure and p is the energy density), corresponding to an ultra-relativistic gas. However, 
the CERN and FNAL colliders provide proton-antiproton collisions with the total 
center-of-mass energy of the order of 1 TeV, corresponding to initial temperatures 
T ~ 10 GeV for which the interaction of the hadrons is strong and has mainly a 
resonant character, the masses of the resonances being comparable with the tem- 
perature. Thus, hadronic matter under these conditions is neither an ideal nor an 
ultrarelativistic gas. 

Corrections in the equation of state due to the interaction of the hadrons have been 
discussed in the literature for the last three decades [1,5,8-16]. These considerations 
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were based mostly on a phenomenological model in which the Landau theory is applied 
to all particles except the leading ones, i.e., the fragments of the initial particles, 
whose characteristics are taken directly from experimental data. The framework for 
the latter considerations is the QCD phase transition of hadronic matter into the 
quark-gluon plasma. 

Hot hadronic matter, in which volumes per particle are a few cubic fermis, is 
certainly made out of individual hadrons. It is clear that, at low temperatures T << 
m w , one has a very rare (and therefore ideal) gas of the lightest hadrons, the pions. 
As the temperature is raised and the gas becomes more dense, one should take into 
account interactions among the particles. This was done using the following three 
approaches: (i) the low-temperature expansion, (ii) the resonance gas, (iii) the quasi- 
particle gas. 

The first approach for the pion gas is based on the Weinberg theory of pion inter- 
actions [O , which uses the non-linear Lagrangian containing all processes quadratic 
in pion momentum. Its first application to calculation of the thermodynamic param- 
eters of the pion gas was done in ref. [1]. One of the important consequences of this 
work was that, after isospin averaging, all corrections quadratic in momenta cancel 
each other, and corrections proportional to the pion mass (Weinberg ittt scattering 
lengths) are nearly compensated, producing negligible corrections at the 1% level. 
Corrections of the second order in the Weinberg Lagrangian were also estimated in 
[1], and a much more systematic study of the problem including quartic terms in 
the mesonic Lagrangian were made in ref. | 18f . It has been shown that taking into 



account the pion rescattering leads to corrections quartic in momenta which, after 
isospin averaging, give a nonzero correction to the pressure, providing a scale for the 
temperatures at which deviations from the ideal pion gas formula become noticeable: 



p = — T 
F 30 



T 
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- int 



150 MeV. 



Without going into discussion of these works we remark that the applicability of this 
approach is limited by temperatures T < 100 MeV, for which typical collision energies 
are significantly below resonance. However, such T are lower than even the lowest 
temperature available in experiments, because the so-called break-up temperatures 
are typically T ~ 120 — 150 MeV. 

The idea of resonance gas was first suggested by Belenky and Landau [plj, as early 
as 1956. They used the Beth-Uhlenbeck method IfZOfl , well known in the theory of the 
non-ideal gases [21 1 , the main idea of which consists in the calculation of the number 



of states within the normalization volumewith the scattering phase shifts taken into 
account. At large distances r between the particles, their wave function is given by 
if>e( r ) ~ sinQor + 8t)/r, and the boundary condition i/jt(R) = at the boundary of 
the normalization volume picks out the states with momenta pR + 5f{p) = ixn. The 
replacement of the summation over n in the statistical sum by an integral over p gives 
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rise to the following expression (Z is only a part of the statistical sum connected with 
the relative motion of the particles): 



1 r°° ( 

Z = -Y,(2i + 1) / dp [R + 

7T £ JO \ 



dp 



e -E[ P )/T_ (2) 



If there is a resonance in the scattering at a certain momentum p , then 5g(p) changes 
rapidly in the neighborhood of Pq. It is easily seen that the contribution of the reso- 
nance to the integral (2) is the term exp (—E(p )/T), the same as for a bound state 
with that energy. Therefore, resonances and stable particles should be taken into 
account on an equal footing in the thermodynamic characteristics of hot hadronic 
matter; e.g., the formulas for the pressure and energy density in a resonance gas read 
(we neglect the chemical potential for simplicity) 

, = £* = £«^E(±U~*^. (3) 

P = Y^Pi, Pi=T < ^-p l , (4) 

where +1 (—1) corresponds to the Bose-Einstein (Fermi-Dirac) statistics, and g>j are 
the corresponding degeneracies (J and / are spin and isospin, respectively), 

(2Jj + l)(2/j + 1) for non — strange mesons 
4(2 Ji + 1) for strange (K) mesons 

2(2J i + l)(2/ i + l) forbaryons 

These expressions may be rewritten with the help of a resonance spectrum, 

p= dmT{m)p(m), p{m) = — — 2^(±1) + = , (5) 

J mi Zty r r 

dm r(m)p(m), p(m) = T ^ - —p(m), (6) 

mi dl 

normalized as ma 

/ dmr(m) =Y,9i, (7) 

J mi i 

where mi and m 2 are the masses of the lightest and heaviest species, respectively, 
entering the formulas (3), (4). 



In both the statistical bootstrap model 0, [K]] and the dual resonance model |l I 
a resonance spectrum takes on the form 



m) ~ m e 



O f, m / T 
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where a and To are constants. The treatment of hadronic resonance gas by means of 
the spectrum (8) leads to a singularity in the thermodynamic functions at T = To 



|S], |TU| and, in particular, to an infinite number of the effective degrees of freedom in 
the hadron phase, thus hindering a transition to the quark-gluon phase. Moreover, 
as shown by Fowler and Weiner |12|], an exponential mass spectrum of the form (8) 



is incompatible with the existence of the quark-gluon phase: in order that a phase 
transition from the hadron phase to the quark-gluon phase be possible, the hadronic 
spectrum cannot grow with m faster than a power. 

In 1972, Shuryak ||, in place of (8), used the simple power parametrization 

r(m) ~ m fc , (9) 

which "describes experiment in the mass region 0.2 — 1.5 GeV for k ~ 3." Since r(m) 
is a rapidly growing function, the main contribution to integrals of type (10) is given 
by the mass region in which exp (m/T) >> 1, and the difference in properties of 
bosons and fermions is irrelevant. Upon substitution of (9) into (5) and (6) one finds 
that the pressure and energy density are proportional to T h+5 . Therefore, the velocity 
of sound c s , defined by 

2 d P nm 

c * = T P > (10) 

turns out to be a temperature independent constant, and the equation of state belongs 
to the class considered by Milekhin ||22|| . For this case, the expressions for energy and 



pressure are as follows, 

o = \T k+5 v = X T k+5 c 2 = 1 (11) 
9 ' P k + A ' s A; + 4' 1 j 

where A is a certain constant. Consequently, for k « 3 ||, c 2 s ~ 0.14. 

In 1975, Zhirov and Shuryak jjTjJ have calculated the velocity of sound, Eq. (10), 
directly from Eqs. (5), (6), and found that at first it increases with T very quickly 
and then saturates at the value of c 2 s ~ 1/3 if only the pions are taken into account, 
and c 2 ~ 1/5 if resonances up to M ~ 1.7 GeV are included, hence favoring k « 1 
in (11) and the equation of state p, p ~ T 6 , p = p/5. Such an equation of state was 
suggested by Shuryak in his book |J of 1988 (and referred to also in [1]), by fitting 
the real resonance spectrum, 

p ~ (20 GeV~ 2 ) T 6 , (12) 



and called the "realistic" equation of state of hot hadronic matter. It gives c 2 s = 0.20, 
in agreement with the theoretical models considered in earlier works |^3|, ^4], pill . 
In ref. ]14| Shuryak, in order to describe the behavior of the energy density with 



temperature p oc T 6 , proposed to consider hot hadronic matter as a gas of quasi- 
particles, which have quantum numbers of the original mesons, but with dispersion 
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relations modified by the interaction with matter (similarly to Landau's idea of "ro- 
tons" , to explain the growth of the energy density of liquid 4 if e with temperature 
more rapidly than T 4 ). Since in such an approach hadrons are included as physical 
degrees of freedom of hot hadronic matter, they are assumed not to be absorbed too 
strongly, i.e., to be good quasiparticles. In ref. |L4| an attempt was made to guess 
what these dispersion relations should be, in order to obtain the expected behav- 



ior of the energy density with temperature. In [16[] these dispersion relations were 
calculated in explicit form. 

It may seem strange that the different approaches (i.e., (i) and (ii) out of the 
three mentioned above, the approach (hi) was invented in order to justify the results 
given by (ii)) should lead to different results, viz., Eqs. (1) and (12); moreover, it 
is not clear which result, if it does, corresponds to the genuine thermodynamics of 
the hadronic resonance gas. In this article we shall show that with the correct form 
of the hadronic resonance spectrum which has been established by the authors in a 
series of papers [26|, |2^, 53], the equation of state of hot hadronic matter takes on the 
form (1), but reduces to (12) in the case of only lower mass resonances being taken 
into account. 



Hadronic resonance spectrum and the equation 
of state 



In our previous work [^6| we considered a model for a transition from a phase of 
strongly interacting hadron constituents, described by a manifestly covariant rel- 
ativistic statistical mechanics which turned out to be a reliable framework in the 
description of realistic physical systems [ffiifl , to the hadron phase described by a 
resonance spectrum, Eqs. (5), (6). An example of such a transition may be a rela- 
tivistic high temperature Bose-Einstein condensation studied by the authors in ref. 
1H, which corresponds, in the way suggested by Haber and Weldon pi], to sponta- 



neous flavor symmetry breakdown, SU(3)f — > SU(2)j x [7(l)y, upon which hadronic 
multiplets are formed, with the masses obeying the Gell-Mann-Okubo formulas |3"2" 



m 



a + bY + c 



T 



(13) 



here / and Y are the isospin and hypercharge, respectively, I is 2 for mesons and 1 
bor baryons, and a, b, c are independent of / and Y but, in general, depend on (p, q), 
where (p, q) is any irreducible representation of SU(3). Then the only assumption on 
the overall degeneracy being conserved during the transition leads to the unique form 
of a resonance spectrum in the hadron phase: 



r(m) = Cm, C = const, 



(14) 
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in agreement with the results of Zhirov and Shuryak ]T3|] found on a phenomenological 
ground. We have checked the coincidence of the results given by a linear spectrum 
(14) with those obtained directly from Eq. (3) for the actual hadronic species with the 
corresponding degeneracies, for all well-established hadronic multiplets, both mesonic 



and baryonic, and found it excellent [26|. Therefore, the fact established theoretically 



that a linear spectrum is the actual spectrum in the description of individual hadronic 



multiplets, finds its experimental confirmation as well. In our recent papers |27|, ^] 
we have shown that a linear spectrum of an individual meson nonet is consistent with 
the Gell-Mann-Okubo mass formula 

ml + 3mg = 4m?/2 (15) 



(in fact, this formula may be derived with the help of a linear spectrum [|27|]), and 
leads to an extra relation for the masses of the isoscalar states, m ' and m », (of 
which belongs to a mostly octet), 

mjy + mjy/ = Wq + mg = 2mw 2 , (16) 

with mi, mi/2, m§, m being the masses of the isovector, isospinor, and isoscalar octet 
and singlet states, respectively, which for an almost ideally mixed nonet reduces to 



27, 28 



m Q n ~ m\, mjy ~ 2m^ 2 — m\. (17) 



We have checked the relation (17) in ref. |28] and shown to hold with an accuracy of 



up to ~3% for all well-established nonets. In ref. |27[] we have generalized a linear 
spectrum to the case of four quark flavors and derived the corresponding Gell-Mann- 
Okubo mass formula for an SU(4) meson hexadecuplet, in a good agreement with 



the experimentally established masses of the charmed mesons. In ref. |33| we have 
applied a linear spectrum to the problem of establishing the correct qq assignment 
for the problematic meson nonets, like the scalar, axial-vector and tensor ones, and 
separating out non-gg mesons. In this paper we shall apply a resonance spectrum to 
the derivation of the thermodynamic characteristics of hot hadronic matter. 

Let us first consider the case of the meson resonances alone. For this case, the 
normalization constant C of a linear spectrum (14) was established in ref. for 



a nonet, one has 9 isospin degrees of freedom lying in the interval (m » ~ mi, m ' 
Therefore, Eq. (7) gives 

f m o' 

C I dm m = 9, 



and hence 



C= , 18 2 ^^27GeV- 
mjy — m^, A 



where the difference A = mjy — is determined by a distance between the parallel 
Regge trajectories for the uj and </> resonances, which are described by the straight 
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lines J = 0.59 + 0.84M 2 and J = 0.04 + 0.84M 2 , respectively, so that A ~ 0.65 GeV 2 . 
We note further that the meson nonets may be arranged in the pairs of nonets which 
have equal parity but different spins (which differ by 2), e.g.|] 

1 3 Pq jpc = 0++; ao (i320), / (1300), / (1525), #5(1430), 

1 3p 2 jpc = 2++) fl2 (i320), / 2 (1270), / 2 (1525), # 2 *(1430), 

1 i D 1 J PC = 1— , p(1700), w(1600), If* (1680), (no <p candidate), 
1 3^ jpc = 3 — ; p3 (i700), u 3 (1600), 03(1850), # 3 *(1780), 

3 % J PC = 0"+, tt(1770), 77(1760), if (1830), (no 77' candidate), 

1 1 D 2 J PC =2-+, vr 2 (1670), if 2 (1770), (no 77,77' candidates), 
and occupy the mass interval of an individual nonet but have 18 isospin degrees of 
freedom in this interval, i.e., twice as much as that for an individual nonet. Moreover, 
as the temperature gets closer to the critical one of chiral symmetry restoration, 
we expect the chiral partners (the states with equal isospin but different parity) 
have equal masses and form parity doublets. The work of DeTar and Kogut [j34fl 
shows convincingly that the "screening masses" of chiral partners are different below 
and become equal above a common T c . This work was carried out for four quark 



flavors. Similar results were obtained for two flavors by Gottlieb et al. [35]. In these 
calculations, the chiral partners were (vr,cr), (p, a\) and (iV(|+), N(^— )). Thus, we 
expect the correct density of states per unit mass interval to be twice as much as that 
for an individual nonet, and hence, the correct normalization constant is 

C ~ 54 GeV" 2 . (19) 

Once the mass spectrum of a nonet (with a given fixed spin) is established to be 
linear, one may take into account different nonets with different spins in Eqs. (3), (4). 
As shown in ref. [p6| , since the particle spin is related to its mass, J« ~ a mf, a 
being a universal Regge slope, the spin degeneracy turns out to be proportional to 
the mass squared, and the account for different nonets results in the following mass 
spectrum, 

r'(m) = C'm 3 , C' = 2a C ~ 90 GeV" 4 , (20) 

which is the actual resonance spectrum of hadronic matter and should lead, through 
(11), to the equation of state 

p,p~T 8 , p = p/7. (21) 



Bebie et al. []36[ have calculated the ratio p/p directly from Eqs. (3), (4), with all 
known hadron resonances with the masses up to 2 GeV taken into account, and 
found that the curve p/p first decreases very quickly and then saturates at the value 
of p/p ~ 7, as read off from Fig. 1 of ref. ]j36[| , in agreement with (21). 

^or the scalar meson nonet, we use the qq assignment suggested by the authors in ref. pq, p3j. 
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In order to show that the obtained normalization constant is correct, we note that 
the number of states with the masses up to M, given by the mass spectrum (20), is 



N(M) = ~ 22.5 (M, GeV) 4 . (22) 

For, e.g., M = 1.25 Eq. (22) gives 

N(1.25) ~ 55. (23) 

The masses up to 1.25 GeV have the members of the pseudoscalar and vector meson 
nonets, and the /ii (1170), &i(1235) and ai(1260) mesons, the mass of the latter was 
indicated by the recent Particle Data Group as 1.23 GeV |57fl. We do not include 



the scalar mesons ao(980) and /o(980) which seem to be non-gg objects ||33|| , but may 
include the /o(1300) meson which has the mass lying in the interval 1 — 1.5 GeV, 
according to the recent Particle Data Group. Thus, we have 9+27+1+9+9+1=56 
actual mesonic species having the masses up to 1.25 GeV, in excellent agreement with 
(23). 

For M — 1.7 GeV, Eq. (22) gives 

N(1.7) ~ 188. (24) 
As seen in the Meson Summary Table the masses up to 1.7 GeV have the 



members of the following nonets: 1 %, 1 3 Si, 1 1 P 1 , 1 3 P , 1 3 Pi, 1 3 ^2, 2 1 S , 2 3 Si. 
Therefore, one has 

(20 spin states) x (9 isospin states) = 180 states, (25) 

in good agreement with the result (24) given by a cubic spectrum. 
For M = 2 GeV, Eq. (22) gives 

N(2) ~ 360. (26) 

The masses up to 2 GeV have the members of all the nonets indicated in [R3 except 



for the 1 3 _F 4 and 2 3 P 2 nonets. In this case, one has 

(41 spin states) x (9 isospin states) = 369 states, (27) 

again in good agreement with the result (26) given by a cubic spectrum. Thus, 
we consider the cubic spectrum (20) as granted by the actual experimental meson 
spectrum. 

Now, as the actual meson resonance spectrum is established, one uses this spec- 
trum in Eqs. (5), (6), in order to obtain the equation of state of hot hadronic 
(mesonic) matter. First we note that one may approximate the particle statistics 
by the Maxwell-Boltzmann one, because of the richness of the spectrum r(m) ~ m 3 
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at large m; then the sum over r in Eqs. (5), (6) may be approximated by 7r 4 /90 = 
C(4) = J2r l/ 7 " 4 ) which is the asymptotic form of this sum for T >> m. One obtains, 
therefore, 

,2 C \2 „ / m y2 



P 



7T 

30 



T 4 + 



180 



T / dm m" 



mV 
TJ 



(2* 



P, 



where we have separated out the contribution of the pions which may well be treated 
as massless at temperatures ~ 150 MeV, and taken into account the remaining parti- 
cle species by an integration over m with the resonance spectrum (20), which therefore 
has the lower limit ~ 0.5 MeV ~ uik- We note that this treatment of the hadronic 
resonance gas corresponds, in view of (3), (5), to a collection of free (non- interacting) 
particles, which is completely justified at lower temperatures, since chiral symmetry 
suppresses the interactions of low energy Goldstone bosons both among themselves 
and with massive hadrons, but becomes an approximation as the temperature in- 
creases. Since the main contribution to integrals of the type (28) is given by the mass 
region in which m » T, one may extend the upper limit of integration to infinity 



and neglect the lower limit, and obtain, through the formula |39 



dx x^K v (ax) = 2"- 1 a^- 1 T 



P 



~ — T 4 
30 



1 + 



1+fl + V 



T 



r 



l + n 



160 MeV 



~ — T 4 
10 



3 V160 MeW 



(29) 
(30) 



The formulas (29), (30) represent the equation of state of hot mesonic matter. 

If one restricts himself to the lower mass resonances (i.e., the lower spin, J = 0, 1, 
nonets) alone, one may use the linear spectrum (14) with C defined in (19) and obtain 



P 



Ctt 2 
180 



dm m 



23 GeV -2 T 



(31) 



in apparent agreement with the Shuryak's "realistic" equation of state (12). 

The formulas (29), (30) have been obtained for the meson resonances alone. There 
is no difficulty of principle to consider the baryon resonances in a similar way, with 
the inclusion of two chemical potentials, for both conserved net baryon number and 
strangeness. As we have checked in ref. ||26|| , a mass spectrum of the SU(3) baryon 



multiplets is linear, as well as for the meson nonets, although to establish its corre- 
spondence to the Gell-Mann-Okubo formulas is more difficult than for a meson nonet, 
since these formulas are linear in mass for baryons (more detailed discussion is given 
in ref. ||26|| ). Recent result of Kutasov and Seiberg |4(| shows that the numbers of 
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bosonic and fermionic states in a non-supersymmetric tachyon-free string theory must 
approach each other as increasingly massive states are included. The experimental 
hadronic mass spectrum shows that in the mass range ~ 1.2 — 1.7 GeV, the number 
of baryon states nearly keeps pace with that of meson states [[|l|] (and, therefore, is 
well described by the same cubic spectrum as for the mesons, Eq. (20)). Above ~ 1.7 
GeV, the number of the observed baryons begins to outstrip that of the mesons, and 
then greatly surpasses the latter at higher energies (indicating, therefore, that the 
baryon resonance spectrum grows faster than (20) in this mass region, since the cubic 
spectrum (20) describes the meson resonances well, up to, at least, 2 GeV, as we have 
seen in Eqs. (23)- (27)). The explanation of this behavior of the experimental reso- 
nance spectrum was found by Cudell and Dienes in a naive hadron-scale string picture 
[4"2"|: the ratio of the numbers of the baryon and meson states should, in fact, oscillate 
around unity, with the mesons favored first, then baryons, then mesons again, etc. 
Keeping in mind this picture, we may assume that the "in-average" baryon resonance 
spectrum has the same form, Eq. (20), as the meson resonance one. If one now 
neglects, for simplicity, the baryon number and strangeness chemical potentials (i.e., 
considers the case of both zero net baryon number and strangeness), and takes into 
account the baryon resonances along with the meson ones by the mass spectrum (20), 
one will obtain the same formulas (28) but with an extra term in the r.h.s. which 
differs from the second term of Eq. (28) by the factor 7/8 (since for the baryons, the 
sum over r in Eqs. (5), (6) is approximated by 77r 2 /720 = 7/8 £(4) = Z) r (— l) r+1 /r 4 ). 
These formulas will further reduce to the relations 



P 



~ —T 4 
30 



' ' 140 Mev) 



(32) 



P 



~ — T 4 
10 



1 + 



7 



T 



3 V140 MeV 



(33) 



which represent the equation of state of hot hadronic (both mesonic and baryonic) 
matter. 

Comparison of Eqs. (29), (32) with (1) shows that the agreement between the 
both approaches, (i) the low-temperature expansion for the pion interacting gas, and 
(ii) the hadronic resonance gas, is very good, if the latter is treated with the help of 
the realistic mass spectrum, Eq. (20). 

We shall use Eq. (1) and the corresponding formula for the energy density, 



P 



~ —T 4 
10 



1 + 



7 



T 



3 V150 MeV 



(34) 



as the equation of state of hot hadronic matter, in agreement with the results (29), (32) 
provided by the hadronic resonance spectrum, Eq. (20). It is seen in (34) that the 
ratio p/T 4 rapidly grows with T, proportionally to T 4 . This behavior is to be expected, 
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since the value p/T 4 = 7r 2 /10 ~ 1 characteristic of a gas of the pions is small compared 
to the value p/T 4 = 47.57r 2 /30 ~ 15.6 which pertains to the plasma of free quarks and 
gluons (47.5 being the number of degrees of freedom in the plasma with N c = Nj = 3). 
The energy density stored in the excited states reaches the energy density of the pionic 
component at T ~ 120 MeV, as follows from (34), in agreement with the value ~ 125 



MeV obtained by Bebie et al. [^6| through direct calculation with the help of Eq. (4). 
The quark-gluon plasma value p ~ 15.6 T 4 is reached at T ~ 235 MeV. However, the 
energy densities of the hadron and quark-gluon phases need not match at the critical 
point, since the phase transition may liberate latent heat. The only pressure in both 
phases should match at the critical temperature, for which, therefore, a naive estimate 
may be obtained by equating (1) with the plasma's p = 47.57r 2 /90 T 4 : T c ~ 290 
MeV. The phase transition from the hadron to the quark-gluon phase is expected 
to occur at much lower temperature. Indeed, the approximation of the hadronic 
resonance gas by a collection of free particles is justified up to temperatures ~ 150 
MeV. Above this point, the interaction among the constituents of the gas rapidly 
grows with T and inelastic collisions become increasingly important. Moreover, direct 



calculation, with the help of Eqs. (3), (4), made by Bebie et al. |36[] shows that 
the mean distance between two particles of the gas reaches the value d — 1 fm at 
T ~ 190 MeV, indicating that Eqs. (3), (4) pack the hadrons very densely there 
and suggesting, therefore, that color cannot remain confined much beyond this point. 
Thus, the equation of state (1),(34) may be an adequate representation of (3), (4) 
up to temperatures ~ 190 MeV. The phase transition to the quark-gluon plasma 
is expected to occur at the same temperatures, ~ 180 — 190 MeV.0 We conclude, 
therefore, that the equation of state (1),(34) is the realistic equation of state of hot 
hadronic matter in the whole temperature range of the latter: — T c , with T c being the 
critical temperature of the hadron to quark-gluon phase transition. In the following 
section we shall build a model for the thermodynamics of this transition. 



3 Phase transition to the quark-gluon plasma 

In the simplest schematic model the quark-gluon plasma is described by the equation 
of state 

47^2 47^2 , R ,o,x 

p= T — B, p = T + B, 35 

y 90 ^30 K ' 

where B is the bag constant, and 47.5 the effective degeneracy of 8 gluons and 3 

massless quarks. In going from the hadron phase (and its nonperturbative vacuum) 

2 The value of the critical temperature of the chiral symmetry restoration transition, to be asso- 
ciated with the hadron to quark-gluon one, calculated with the contribution of the massive states 
being taken into account, is Q] ~ 190 MeV if only the meson resonances, and ~ 175 MeV if both 
the baryon and meson resonances are considered, in agreement with the results obtained previously 
by Gerber and Leutwyler |lq|. 
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to the deconfined quark-gluon phase (with perturbative vacuum) we restore conformal 



invariance. Consequently [44], we should use the value of B given by the conformal 
anomaly: 



B 



■(0|^|0) = -^(0|^(G^) 2 |0). 



(36) 



Here 9£ is the trace of the energy-momentum tensor and (0|(G^) 2 |0) is the gluon 
condensate. The /3-function of the renormalization group is given by |E3] 



16tt 2 
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(37) 



for Nf flavors. The value of the gluon condensate, 



471-2 



<0|(GL ) 2 |0) ^ (330 MeV) 4 . 



(38) 



is well determined by the states of charmonium using QCD sum rules [JB]. For three 
flavors, we then find 

B l/i ~ 240 MeV, (39) 

and Eqs. (35) reduce to 



V 



47.5tt 2 
90 



rp4 



1 



^160 MeV^j 



1 i 



47.5vr 2 . 



P 



30 



rp4 



1 /160 MeVV 

i — y — ) 



. (40) 



If the hadron phase is described, e.g., by Eqs. (29), (30), a simple equating of the 
pressure in both phases (Eqs. (29) and (40)) gives two real positive roots: 



7\ « 166 MeV, T 2 w 308 MeV. 



(41) 



In the temperature range 7\ < T < T 2 , the curve p = p(T) for the quark-gluon phase 
goes above the corresponding curve for the hadron phase, but at T > T 2 , the situation 
changes: the quark-gluon curve goes below the hadron one. In this picture, therefore, 
the hadron phase is again thermodynamically favored at high temperatures, and the 
existence of the quark-gluon phase is restricted to a limited range of temperatures, 
which is clearly unphysical. The reason for such unphysical behavior predicted by 
Eqs. (29), (40) is, except for the applicability of Eqs. (29), (30) limited by T ~ 190 
MeV, the effective degeneracy in the hadron phase which grows, according to (29), 
as g e ff ~ T 4 . As shown by Brown et al. J44|, the free energy of an interacting system 
of massless microscopic constituents (and hence, of the effective degrees of freedom 
which are intended to mimic these interactions) should be greater than the free en- 
ergy of the noninteracting microscopic constituents at any temperature. When these 
interactions are realized in terms of massless hadrons, this constraint implies that the 
number of the effective (hadronic) degrees of freedom should always be less than 47.5 
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and that it should approach 47.5 in the high-temperature limit. In other words, the 
number of degrees of freedom in the effective excitations (hadrons) should not exceed 
that of the underlying quarks and gluons. Specifically, the hadronic resonances are 
intended to simulate attractive channels of the microscopic constituents. The stability 
of the trivial microscopic vacuum at zero temperature against the production of in- 
finitely many particles requires the existence of repulsive channels for the microscopic 
constituents, which is not taken into account in Eqs. (1),(29),(32). At the one-body 
level, such repulsion is most easily simulated by associating an excluded volume with 
each hadron [44 1 . However, the most naive assumption of a temperature independent 
excluded volume represents an "overkill" in the high-temperature description of the 
now massless hadrons, since it leads to an effective degeneracy which vanishes as the 
temperature grows up [P]. A more realistic description of the hadron phase can 
be obtained using an excluded volume which decreases like 1/T 3 for large T, which 
enforces the limitation of the number of the effective degrees of freedom by a fixed 
number, which may be made taken a desired value (e.g., 47.5) by the corresponding 
choice of the parametrization for a temperature dependent excluded volume ||44[ . An 
alternative approach is to restrict ourselves to a finite number of hadronic states, 
which should be chosen to be 47.5, since, according to the arguments given above, we 
should have 47.5 effective massless degrees of freedom in the hadron phase at high 
temperature. For our present semiquantitative purposes, this approach provides an 
adequate alternative to the use of excluded volumes. 

The most natural way to restrict the effective degeneracy in the hadron phase 
(to 47.5) is truncating the resonance spectrum (20) when the number of the effective 
hadronic degrees of freedom is 47.5. This makes certain sense from the physical point 
of view, similar to the vibrations of a Debye solid where the phonon spectrum is 
truncated when the number of phonons is 3N, with N being the number of atoms. It 
is seen in (22) that M ~ 1.2 GeV corresponds to iV = 47.5; we consider, therefore, the 
following expression for the pressure in the hadron phase with 47.5 effective degrees 
of freedom: 

P = ™ T + ^7T T / dmm 3 ( — ) K 2 ( — ) , (42) 



30 180 Jo.5 \TJ *\T. 

where we have, as previously, separated out the contribution of (massless) pions. It is 
now seen that at T ~ 1.2 GeV, this formula reduces, through K 2 (x) ~ 2/x 2 , x « 1, 
to 

7r 2 C'n 2 r 1 - 2 

p = !L T 4 + hl^ T 4 f dm m 3 (43) 

y 30 90 Jo.b V ; 

since C' J Q f dm m 3 = 44.5 (= 47.5 — 3 of the pions), it follows from (43) that at high 
enough temperature, 

47.5tt 2 . . 

P = T\ 44 

y 90 V ; 

i.e., we have 47.5 effective massless degrees of freedom in the hadron phase. Fur- 
ther comparison of Eqs. (40) and (44) shows that the transition to the deconfined 
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quark-gluon phase is still prohibited by the bag constant B : the pressure curve (40) 
goes below the corresponding (44) at high T , and hence the hadron phase is thermo- 
dynamically favored there, even if the curves (40) and (44) intersect at some lower 
temperature. Thus, one sees that even in the model with the restricted effective 
degeneracy in the hadron phase, the deconfinement transition is absent. This con- 
clusion is in agreement with lattice gauge calculations which show the only phase 



transition, viz., the restoration of chiral symmetry [34, pq], not deconfinement. As 



stressed out in ref. f44]| , the deconfinement transition for pure glue, used in lattice 



calculations as a model problem for many years, is a fictitious transition; it is com- 
pletely changed by the introduction of light mass quarks. Moreover, the results of 
lattice gauge calculations indicate that the gluon condensate is essentially unchanged 
during the chiral restoration transition WJ\ . Even at a temperature as high as ~ 290 



MeV, lattice calculations find about half as much gluon condensate as at zero tem- 
perature [Q. Apparently, the conformal anomaly (and the bag constant associated 
with the former) is relatively unaffected by what happens to the quark condensate. 
The same conclusion was reached by Brown |49| in the connection with the variations 



of a chemical potential rather than temperature. Thus, even at temperatures well 
above the critical one of chiral symmetry restoration, T x (~ 190 MeV), gluons remain 
condensed. Consequently, although the mesons can be regarded as quark- ant iquark 
pairs, each quark must be connected with an antiquark by a "string" (i.e., a line 
integral of the gauge vector potential), in order to preserve gauge invariance. It is 
difficult to include the important consequences of such quark-antiquark correlations 
in the thermodynamics of the gluon condensed system, unless one keeps regarding the 
mesons as the correct effective degrees of freedom even above T x . Since our equation 
of state (42) does regard the mesons as the effective degrees of freedom at any T, we 
may model the chiral symmetry restoration transition with the help of Eq. (42). 



Campbell et al. |50[ have shown how to build the correct operator scaling proper- 
ties into a low-energy effective Lagrangian in the Nambu-Goldstone (chir ally- asym- 
metric) sector. Brown and Rho [[51] have noted that this implies (in the chiral limit 
in which the bare (current) quark masses are zero) that 

.. = $(T), (45) 



m p m N f w 



where $(T) describes the common scaling of all hadronic properties. Here the starred 
quantities denote in-medium values at finite temperature and/or density. The order 
parameter for chiral symmetry breaking is the effective pion decay constant, /* [[)(], 



51~| . One can associate f*^0 with with the melting of the quark condensate, (qq)*, 



with increasing temperature. The form of $(T) was justified in ref. [ 5"T |: 
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We shall, however, use a different form of $(T), viz., 



$(T) 



(9?) ' 



(47) 



suggested by the Nambu model and shown by Koch and Brown |52] to model lat- 
tice gauge data better than (46). With the temperature dependence of the quark 



condensate found by Gerber and Leutwyler |18j and confirmed in our recent paper 

(qq)* 



43 



(??) 

Eq. (45) leads to the relation 



1 



T ~ 
± x — 



190 MeV, 



(4£ 



m 



m 1 



(49) 



for all hadronic species. With (49), Eq. (42) takes on the form 



V 



TV 



30 180 



/■1.2 

T 4 / dm m 3 
Jo.5 



m 



m 



(50) 



which upon the integration reduces, through C'(1.2 4 — 0.5 4 )/4 = 44.5, to 



P 



7T 

90 



rpi 



3 fm^ 
2 V^T 



T 



+ 89 



y 5 K 5 (y) - x 5 K 5 (x) - 67/ 4 K 4 (y) + 6x 4 K 4 (x) 



x 1 



X = 



1.2$(T) 



y 



0.5$(T) 



(51) 



where $(T) is given in (47), and we have recovered the pion contribution term with 
the finite pion mass. Temperature dependence of the effective degeneracy, g e ff = 
p/Psb, with p SB = 47.57r 2 /90 T 4 , as well as those of p/psb, with p S B = ^Psb, 
and p/p (~ c 2 ), are shown in Figs. 1-3, respectively. One sees that in the model 
considered here, chiral symmetry restoration corresponds to a smooth crossover in 
the thermodynamic variables p and p. Let us compare these results with those of 
lattice gauge calculations. To date the most complete lattice gauge calculations on 
energy production for temperatures near T x are those obtained by Kogut et al 



For calculations involving two low mass quarks, they find rapid increases in both 
the quark and gluon energy densities beginning at (3 = 5.30 on a 6 x 12 3 lattice 
(j3 is related to the inverse coupling constant, j3 = 6/g 2 ). By ft = 5.35, the sum of 
the quark and gluon energy densities is essentially blackbody (and corresponds to 
that of an ultrarelativistic Stefan-Boltzmann gas), although in the case of gluons the 
accuracy of their calculation is not good enough. They find "a smooth p^d/T 4 curve 
with considerable suppression [compared to blackbody] in the transition region." The 
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measured transition is a smooth crossover, in contrast to a sharper change obtained 
previously with smaller lattices. Earlier calculations of the behavior of the quark and 
gluon energy densities in this temperature region are summarized by Petersson |3J] 
and, within their accuracy, are roughly consistent with the results of Kogut et al. 
53 . 

This smooth crossover is in agreement with the results obtained in this paper 
for the model of the hadron to quark-gluon transition with the truncated resonance 
spectrum and the temperature dependent effective hadron masses. It is seen in Figs. 
1-3 that the crossover begins at T ~ 140 MeV (corresponding to a dip in the In- 
dependence of p/p at those temperatures seen in Fig. 3 [0]), i.e., at T ~ 3/4 T x , and 
finishes at T ~ T x , in agreement with the predictions made by Brown et al. |[44|| . At 
T > T x ~ 190 MeV, the system is essentially blackbody. 

Although the results obtained in this paper are semiquantitative, they agree well 
with available lattice gauge data and reflect the main features of the latter |] i.e., a 
smooth crossover in the thermodynamic variables taking place in a temperature range 
~ 50 MeV, and the temperature dependent ratio p/psb approaching unity from above 
as the temperature increases. 

4 Concluding remarks 

We have obtained the equation of state of hot hadronic matter, by taking into account 
the contribution of the massive states with the help of the resonance spectrum r(m) ~ 
m 3 justified in our previous papers. This equation of state is in agreement with 
that provided by the low-temperature expansion for the pion infracting gas. We 
have shown that in this picture the deconfinement phase transition is absent, in 
agreement with lattice gauge calculations which show the only phase transition of 
chiral symmetry restoration. The latter was modelled with the help of the restriction 
of the number of the effective degrees of freedom in the hadron phase to that of the 
microscopic degrees of freedom in the quark-gluon phase, through the corresponding 
truncation of the hadronic resonance spectrum, and the decrease of the effective 
hadron masses with temperature, predicted by Brown and Rho. The results are in 
agreement with lattice gauge data, thereby confirming the prediction of Brown and 
Rho, and show a smooth crossover in the thermodynamic variables which begins at 
T ~ 3/4T x and finishes at T ~ T x , where T x ~ 190 MeV is the critical temperature 
for the chiral symmetry restoration transition, i.e., takes place in a temperature range 
~ 50 MeV. At T = T x the system ends up in the chirally symmetric phase and is 
essentially blackbody. 

3 These features are seen, e.g., in Fig. 1 of the most recent review on lattice gauge calculations |5q| : 
the crossover starts at T ~ 150 MeV and finishes at T ~ 200 MeV, the ratio pj psb approaches unity 
from above as the temperature increases. With the account for the error bars of this calculation, 
the system is essentially blackbody at T just above 200 MeV. 
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Although our results are in good agreement with lattice gauge data, we, along 
with the authors of rcf. [11] . have at present no way of how to decompose the hadron 
pressure (and energy density), 

2 

p = 9 -^~T\ g eff - 47.5 with T ^ oo, 

into the separate quark and gluon components both above and below T x . Also, we 
are not clear about what happens above T x in a real world. This is the matter of 
subsequent investigations. 
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FIGURE CAPTIONS 



Fig. 1. The ratio p/psb, with p given in Eq. (51), as a function of temperature. 



Fig. 2. The ratio p/psB as a function of temperature. 



Fig. 3. The ratio p/p as a function of temperature. 
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